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by 

V.I.  Ollker  and  P.  Haltaan 


RESEARCH  PROPOSAL  AND  PROGRESS  REPORT 


J  The  physical  phenomena  described  by  nonlinear  partial  differential 
equations  have  become  at  present  the  central  theme  of  investigations  by  many 
researchers.  A  good  understanding  of  most  physical  processes  requires 
accounting  for  nonlinear  effects  and,  consequently,  methods  for  studying 
nonlinear  equations  have  to  be  developed. 

Among  nonlinear  equations  the  Dlrlchlet  problem  for  the  Monge-Ampere 
equation  is  the  model  case  for  fully  nonlinear  equations.  The  problem  is 
formulated  as  follows. 


In  Euclidean  plane  R  with  Cartesian  coordinates  x,y  consider  a  bounded 
domain  S2,  a  nonnegative  function  f  :  H  10,«>),  and  a  continuous  function 
^  R.  It  is  required  to  investigate  solubility  of  the  problem 


M(z)  ■ 


3y^  \ 


2  Y 


3  z 
dxdy  J 


f  in  R, 


(1) 


<t»,  r  -  an. 


(2) 


The  equation  (1)  is  perhaps  the  most  simple  representative  of  the  class  of 
nonlinear  equations  of  Monge-Ampere  type  (see  ICHj,  p.  324,  iGTj,  ch.  17).  Such 
equations  have  been  studied  by  many  authors  lA|,  iBj,  [CNSj,  tCYj,  lN|,  [Oj, 
[Pl],  [ P2 ] ,  mainly  in  connection  with  problems  of  existence  and  uniqueness  of 
surfaces  with  prescribed  metric  or  curvature  functions.  However,  they  also  have 
other  important  applications.  In  particular,  the  leading  term  in  the  "balance 


M-r 


Cm 


. i 
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equation"  in  dynamic  meteorology  has  the  form  (1)  iH].  In  a  more  complicated 
form,  an  equation  of  this  type  appears  in  the  von  Rarman  system  of  equations  for 
elasticity  and  also  in  an  inverse  problem  of  geometric  optics  [W];  see  below 
part  B.  It  also  turned  out  that  recent  progress  in  the  study  of  fully  nonlinear 
equations  became  possible  after  Important  properties  of  the  equation  (1)  were 
discovered;  [GTj,  ch.  17.  For  quaailinear  elliptic  and  parabolic  equations  the 
use  of  equations  of  the  form  (1)  is  crucial  in  obtaining  C°  and  Holder 
estimates  (see  [GT],  ch.  17,  [KJ,  and  further  references  there). 

In  spite  of  the  increasing  number  of  papers  in  this  area  the  theory  of  the 
problem  (1),  (2)  and  its  generalizations  Is  far  from  being  complete. 

Researchers  in  the  USA,  USSR,  Germany,  England  and  other  countries  at  present 
are  actively  pursuing  this  direction. 

The  proposers  were  fortunate  to  start  their  research  at  a  relatively  early 
stage  of  all  these  developments.  The  current  project  funded  by  the  AFOSR  Grant 
84-0285  covers  two  major  areas: 

A.  Investigation  of  numerical  methods  for  solving  problem  (I),  (2)  and  its 
generalizations . 

B.  Investigation  of  solubility  of  a  Monge-Aropere  equation  arising  in 
shaped  antenna  design. 

We  describe  now  briefly  the  results  obtained  so  far,  current  activities, 
and  proposed  directions  for  future  studies. 

Part  A.  Huaerlcal  Methods.  The  type  of  the  operator  M  depends  on  the 
function  on  which  it  is  evaluated.  For  that  reason,  one  usually  seeks  a 
solution  of  (1),  (2)  in  the  class  of  functions  on  which  M  has  a  fixed  type,  for 
example,  elliptic  or  hyperbolic.  In  the  case  under  consideration,  the 
requirement  f  >  0  in  Ji  forces  any  function  satisfying  (1)  to  be  an  "elliptic" 
solution,  that  is,  in  the  class  where  M  is  elliptic. 


3 


In  view  of  the  Important  practical  applications  several  heuristic 
approaches  were  suggested  for  numerical  solution  of  some  modified  forms  of  (1), 
(2),  [Ar],  [Sh],  IShlJ.  Though  no  rigorous  analysis  of  these  methods  exists, 
one  may  note  that  they  all  are  local  methods  based  on  a  finite  difference 
approximation  and  linearization.  Because  of  the  strong  nonlinearity  of  N  this 
approach  might  be  successful  only  In  a  neighborhood  of  the  true  solution  and 
therefore,  If  a  priori  a  good  Initial  approximation  Is  not  available,  these 
methods  will  not  produce,  generally,  a  sequence  converging  to  the  true  solution. 

We  Investigated  this  problem  In  detail  and  obtained  the  following  results; 

•  A  special  discretization  scheme  for  (1),  (2)  was  suggested  different 
from  standard  finite  element  or  finite  difference  schemes.  It  can  be 
shown  that  latter  ones  In  known  forma  will  not  work  here; 

•  an  iterative  method  has  been  developed  for  solving  the  discretized 
version  of  (1),  (2); 

*  the  question  of  finding  an  Initial  approximation  In  our  scheme  Is 
completely  and  effectively  resolved;  It  Is  just  a  routine  step  of 
the  Iteration  process; 

*  the  Iterations  are  selfcorrectlng; 

*  global  convergence  la  established; 

’  our  algorithm  is  suitable  for  a  parallel  computer; 

•  a  computer  code  has  been  written  and  tested. 
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The  experience  gained  In  testing  our  procedure  on  different  types  of 


examples,  Including  ones  with  large  gradients,  shows  that  Its  most  effective  use 


will  be  In  combination  with  some  fast  Newton-type  scheme.  More  precisely,  a 


particular  criteria  exists  for  checking  when  the  current  approximation  can  be 


used  as  the  beginning  step  for  a  converging  Newton-type  Iterative  procedure.  In 


this  combined  scheme  we  proved  convergence  (quadratic)  of  the  Newton  Iterates, 


but  the  above  mentioned  criteria  Involves  some  heuristic  arguments  and  more  work 


needs  to  be  done  here. 


The  computer  code  for  the  method  Is  quite  sophisticated;  It  Involves,  as  a 


step,  construction  of  a  convex  hull  of  seta  of  points  In  R  .  There  are  here 


different  approaches  and  the  effectiveness  of  the  algorithm  depends  on  It.  We 


have  been  currently  testing  various  schemes  and  the  code  presently  Is  a 


substantial  Improvement  over  Its  original  version  of  1984-85.  A  full  time 


employed  graduate  student  In  Computer  Science  with  a  Ph.D.  in  Mathematics  has 


been  helping  us  with  this  part  of  work. 


The  results  have  been  submitted  for  publication  (see  the  subsection  on 


publications  and  presentations). 


In  the  third  year  of  the  funding  period  we  Intend  to  study  the  already 


mentioned  before  "combined  schemes",  two  sided  approximations,  and  extensions  to 


equations  of  the  form 


a^z  a^z  f  ^  +  2h^  ) 


/  %  .  9z  dz.  ^ 

+  c(x,y)  — 7  ♦  g(x,y,z,— ,  — )  -  0. 


We  have  already  initiated  work  in  these  directions  and  are  quite  sure  of 
its  sucessful  coropletioi  tie  last  equation  includes  as  a  special  case  the 
balance  equation  of  the  mospherlc  dynamics.  If  time  and  resources  permit  we 
want  to  tailor  the  code  we  have  for  a  parallel  machine.  At  the  University  of 
Georgia  in  Athens  a  Cyber  205  is  available  and  we  have  made  some  contacts  with 
the  program.  We  have  also  established  contact  with  Professor  D.  Hoffman  from 
the  University  of  Massachusetts  at  Amherst  who  agreed  in  principle  to  help  us  to 
analyze  our  work  with  computer  graphics.  His  experience  in  this  area  is  well 
known  [HJ  and  we  hope,  it  will  be  valuable  for  us. 

Part  B.  Applications  to  Shaped  Antenna  Dealgn.  In  a  practical  problem  of 
shaped  antenna  design  it  is  required  to  determine  a  reflecting  surface  such  that 
for  a  given  point-source  of  light  the  reflected  rays  cover  a  prescribed  region 
of  the  far  sphere  and  the  density  of  the  distribution  of  reflected  rays  is  a 
prescribed  in  advance  function  of  the  direction.  It  is  assumed  that  the  power 
density  of  the  source  as  well  as  the  apperture  of  the  Incident  ray  cone  are 
known,  and  the  reflection  process  obeys  the  laws  of  geometric  optics  (see  Figure 
1  on  next  page) .  In  this  form  the  problem  was  posed  by  Westcott  and  Norris 
IWNj  and  later  it  has  been  considered  in  Brickell,  Harder,  and  Westcott  [BMWj, 
Norris  and  Wescott  [NW],  Wescott  [Wj.  Further  references  can  be  found  in  [W], 
The  research  of  these  investigators  has  been  supported  by  Plessey  Radar,  Ltd. 
for  many  years. 

The  problem  admits  a  precise  mathematical  formulation  and  in  this  form  it 
reduces  to  solving  the  equation 


{|VP|^  +  P^)^  det(e^^) 


p  in  w. 


(  3) 


nonlinear  boundary  condition  on  9u, 


(4) 


OBJECT:  To  Illuminate  Target  with  Prescribed  Intensity 
(Through  Given  Aperture) 


Equivalent  to  Determining  the  Surface  and  Position  of  Reflector 
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with  respect  to 
problem;  here  (e 


the  unknown  function  p( >0)  naturally  associated  with  the 
^ j )  is  the  matrix  of  the  first  fundamental  form  e  of  the  unit 


sphere  S  ,  V  the  gradient  in  the  metric  e,  7 


ij 


second  covariant  derivatives 


in  e.  The  condition  (4)  is  somewhat  complicated  to  be  presented  here  without 
considerable  expansion.  In  [NW],  and  also  in  [W],  p.  40,  the  authors  state  that 
the  question  of  existence  of  solutions  is  open.  Conditions  for  uniqueness  were 
given  by  Harder  [M].  In  [NW],  [W]  results  of  numerical  studies  and  applications 
are  presented.  Still  there  are  no  rigorous  convergence  results,  and  even  the 
linearized  version  of  the  problem  has  not  been  investigated.  In  1957  J.  Keller 
[kJ  obtained  some  results  pertaining  to  the  radially  symmetric  case  without 
satisfying  a  particular  boundary  condition. 

With  the  support  of  AFOSR  we  started  our  Investigation  of  the  problem  with 
the  radially  symmetric  (r.s.)  case,  that  is,  when  the  incidence  ray  cone  12  and 
the  far  field  domain  u  are  circular,  the  prescribed  density  of  reflected  rays  is 
a  function  of  the  aslmutal  angle  only,  and  the  reflecting  surface  is  sought  as  a 
surface  of  revolution.  The  following  results  have  been  obtained; 

•  The  problem  splits  naturally  into  two  parts.  In  the  first  part  one 
finds  a  class  of  surfaces  for  which  the  reflected  directions  cover 

the  prescribed  domain  u.  Those  surfaces  can  be  conveniently 


parametrized  by  points  of  u.  In  the  second  part  one  seeks  in  the  above 
class  a  particular  surface  for  which  the  density  of  the  reflected 
rays  is  a  prescribed  function  in  to,  and  which  projects  onto  the  given 

domain  {2; 

2  - 

•  it  la  shown  that  for  any  function  p  £  C  (lo),  p  >  0  in  w,  the  surface  F 


defined  by  Che  map 
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-  r  -  Vp  +  (p  -  p)y,  yew, 


p  -  (p^  +  |Vp|^)/2p  -  |r|.|7p|^  -  <Vp,Vp>, 


satisflea  the  requirement  of  being  a  reflector; 


simple  and  verifiable  necessary  ar.d  sufficient  conditions  for 


solubility  of  (3),  (4).  It  was  shown  that  radially  symmetric 


solutions  can  be  constructed  explicibly  whenever  all  of  the  parameters 


are  appropriate. 


In  a  separate  study  that  Just  has  been  completed  we  proved  that 


nonradially  symmetric  solutions  to  the  problem  for  general  distribution 


densities  which  are  close  (in  C  ’  -  norm)  to  r.s.  solutions  can 


be  found  under  the  natural  restriction  corresponding  to  the 


conservation  law.  The  result  appears  to  be  quite  delicate,  because. 


in  general,  the  problem  is  unstable  and  even  a  small  perturbation  of 


the  data  may  produce  a  situation  where  a  solution  does  not  exist. 


During  the  past  year  we  have  also  experimented  with  numerical  methods  for 


solving  (3),  (4)  for  ^  sufficiently  close  to  a  r.s.  symmetric  density.  In 


r.s.  case  the  equation  is  singular  at  the  endpoint  corresponding  to  axis  of 


revolution,  and  as  a  result  straight  forward  linearizations  about  r.s.  solutions 


do  not  seem  to  work  very  well,  though  in  iWj  some  success  was  claimed. 


During  the  third  year  we  Intend  to  continue  investigation  of  the 


following  topics: 


Solubility  of  (3),  (4)  with  the  right  hand  side  subject  only  to  the 


conservation  of  energy  requirement  (a  necessary  condition  in  this 


problem) ; 


^  ^  ^  'j*  ^ . 


•  numerical  solution  of  (3),  (4)  for  ^  close  to  r.s.  data.  With  our 


I 


latest  work  now  complete  we  have  much  better  Insights  into  the  problem 
and  anticipate  definite  progress  here. 

•  Application  of  numerical  methods  developed  in  Part  A  to  the  problem  (3), 
(4). 

•  Extension  of  the  techniques  used  in  the  solution  of  the  r.s.  case  to 
more  general,  nonlinear  eigenvalue  problem  for  ordinary  differential 
equations . 


: 
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Publlcatlona  and  Preaentatlona 

The  results  obtained  in  Part  A  were  presented  in 

■  a  paper  "On  the  numerical  solution  of  the  equation 

ill  ill  ( ili-\ 

3x2  3y2  3x3y  J 

and  its  discretizations,  I".  It  has  been  submitted  to  the  SIAM 
Journal  of  Numerical  Analysis  and  is  currently  under  revision; 

*  a  lecture  at  a  special  session  of  the  Southeastern-Atlantlc  Regional 
Conference  on  Differential  Equation,  Atlanta,  October,  1985 

*  a  plenary  lecture  was  delivered  on  the  subject  by  one  of  the  principal 
investigators  at  the  South  Eastern  Regional  SIAN  Conference,  March, 
1986. 

The  results  in  Part  B  were  presented  in 

*  a  paper  "Radially  Symmetric  Solutions  of  a  Nonge-Ampere  Equation 
arising  in  a  reflector  mapping  problem",  pp  1-19;  accepted  for 
publication  in  Proceedings  of  UAB  Conference  on  differential 
equations  and  mathematical  physics.  Springer; 

’  a  paper  "Near  radially  symmeytric  solution  of  an  inverse  problem  in 
geometric  optics",  pp  1-23,  to  be  submitted; 

‘  a  report  "On  the  Monge-Ampere  equation  arising  In  the  reflector  mapping 
problem".  Institute  for  Mathematics  and  Its  Applications,  University 
of  Minnesota,  preprint  series  <^198,  pp.  1-43  (this  Is  an  expanded 
version  of  the  first  paper  mentioned  In  this  part); 
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|i  •  an  hour  lecture  given  at  the  Technical  Unlveraity  In  Berlin,  West 

*  Gernany; 

i: 

^  *  a  30->nilnute  presentation  at  the  conference  on  Differential  Equations  and 

^  Mathematical  Physics,  University  of  Alabama,  Birmingham. 


